First-principles structural optimization and electronic structure of 
the superconductor picene for various potassium doping levels 
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We theoretically explore the crystal structures of picene, for which a new aromatic supercon- 
ductivity has recently been discovered for x = 3, by systematically performing first-principles full 
structural optimization covering the concentration range x — 1-4. The crystal symmetry (space 
group) of the pristine picene is shown to be preserved in all the optimized structures despite sig- 
nificant deformations of each picene molecule and vast rearrangements of herringbone array of 
molecules. For K^cpicene (x = 1-4) optimization indicates that (i) multiple structures exist in some 
cases, and (ii) dopants can enter in the intralayer region as well as in the interlayer region between 
the stack of herringbone structures. In the electronic structure obtained with the local density 
approximation for the optimized structures, the rigid-band approximation is invalidated for multi- 
fold reasons: the dopants affect the electronic properties not only through the rearrangement and 
distortion of molecules, but also through hybridizations between the molecules and metal atoms. 
As a consequence the resultant Fermi surface exhibits a variety of multiband structures which take 
diverse topology for Ki picene and K3 picene. 

PACS numbers: 74.20. Pq, 74.70. Kn, 74.70.Wz 



I. INTRODUCTION 

The discovery of superconductivity in potassium- 
doped solid picene, K^^ picene, with Tc = 7-18 K for a 
doping level x 3, recently reported by Mitsuhashi et 
al. [l|, is seminal as the first aromatic superconductor. 
This is expected to open a new avenue in quest for or- 
ganic superconductors. The occurrence of superconduc- 
tivity in the doped picene should not be an exceptional 
case, since the discovery has been followed by a more re- 
cent report of superconductivity in another of aromatic 
compounds, coronene, by Kubozono et al. Q with Tc 
up to 15 K for coronene with x c:^ and also in 
doped phenanthrene with Tc up to 5 K for x 3 by 
Wang et al. [3]. These aromatic compounds are now 
forming a new class in the family of carbon-based su- 
perconductors that include graphite-intercalation com- 
pounds (GICs) and doped fullerenes. For GIC, KCs, the 
first carbon-based superconductor with Tc = 0.1 K[4], is 
followed by CaCe with the highest Tc = 11.6 K among 
GICs[5[. Doped fullerene K3C60, a superconductor with 
Tc = 20 K 0, is followed by Cs2RbC6o with Tc = 33 
K0, and CssCeo with Zi, = 40 K under 15 kbar[8j or 
Tc = 38 K under 7 kbarQ. 

The electronic structures of pristine and doped picene 
were first obtained by the present authors [10] as an ini- 
tial step to elucidate the mechanism of superconduc- 
tivity. It has been revealed that the conduction band 
group is derived from the lowest unoccupied molecular 
orbital (LUMO) and LUMO+1 of a picene molecule. 
The conduction-band structure of solid picene, where the 
molecules are arrayed in a herringbone structure with 
two molecules per unit cell, is indeed shown to be accu- 



rately reproduced by a four-orbital tight-binding model 
downfolded in terms of the maximally-localized Wannier 
functions [ll|. We have also determined the dopant po- 
sitions for Kapicene, where we have employed the exper- 
imental lattice parameters. Andres et al. [lH, on the 
other hand, provided two, fully-optimized geometries of 
Kapicene, one with a herringbone structure and the other 
with a laminar structure. Giovannetti and Capone[l3| 
performed an electronic structure calculations with dif- 
ferent exchange-correlation functionals, and found that a 
hybrid functional leads to an antiferromagnetic state in 
Kapicene. Kim et al. [ij] demonstrated that an antiferro- 
magnetic state is stabilized in Kapicene when the volume 
is increased. They, as well as Giovannetti and Capone, 
suggested that a strong electron correlation exists in the 
system. This is reasonable in that the electronic correla- 
tion energy in solid picene is estimated to be /7 ^i:^ 0.85 
eV with the momentum-dependent electron-energy-loss 
spectroscopy, [l5| which is actually greater than the width 
of the LUMO- and LUMO+l-derived conduction bands. 
This is also the case for potassium-doped picene, as will 
be demonstrated below. 

Since the report of the superconductivity in potassium- 
doped picene, a few theoretical studies aiming at clari- 
fying its mechanism of superconductivity has been re- 
ported. These include phonon studies such as Kato 
et al. [1^, who have estimated the phonon frequencies 
and the electron-phonon coupling for various hydrocar- 
bon molecules. For picene they found that the electron- 
phonon coupling is large enough to attain Tc 10 K. 
Subedi and Boeri[l3 calculated the electron-phonon in- 
teractions in doped picene within the rigid-band approxi- 
mation, and obtained coupHng constants large enough to 
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explain the experimental Tc ^ 18 K in Kapicene. Casula 
et al. [m performed a more elaborate calculation, and 
obtained the electron-phonon coupling large enough to 
explain the experimental Tc ^ 7 K, to which they found 
the dopant and intermolecular phonon modes contribute 
significantly. 

However, one essential question is: what is so special 
about Kspicene as opposed to other values of xl The 
question is important, since it should be intimately re- 
lated to another question: to what extent the aromatic 
superconductivity is general. We can for instance re- 
member that the problem of why Asfullerene {A: alkali 
metal elements) is special played an important role in 
elucidating its electronic properties. One experimentally 
unresolved situation is that, in the doped picene, both its 
atomic structure and even the number of dopant atoms 
actually accommodated in the samples have yet to be 
experimentally determined, although there are some es- 
timates as discussed in Ref.[2]. Given the circumstance, a 
theoretical prediction for the geometry of doped picene is 
thus desirable, especially for systematically varied dop- 
ing levels. This will also give a starting point to dis- 
cuss whether the superconductivity is driven by a phonon 
mechanism, or by an electronic mechanism as in some 
classes of organic superconductors [19j. 

This is exactly the purpose of the present work. For 
organic materials, not only the atomic positions but also 
the unit cell should be relaxed for a reliable optimization 
of the geometry. Hence we have performed here the full 
structural optimization of Ka^picene, covering the concen- 
tration range x = 1-4. 



II. METHOD 

For the electronic structure calculation and the struc- 
tural optimization in the present work, we adopt the 
density- functional theory with the projector augmented- 
wave (PAW) method ^] with the Quantum MAterials 
Simulator (QMAS) package [2l| in the local density ap- 
proximation (LDA)|22|. The Bloch pseudofunctions are 
expanded in plane waves, with an energy cutoff of 40 
Ry throughout the present work. All the self-consistent 
field calculations are performed with 4 x 4 x 4 /c points. 
The structural optimization allows the atomic coordi- 
nates and the cell parameters to be relaxed with no con- 
straint imposed on them. In each iteration of the self- 
consistent field calculation, the forces and stress [23| felt 
by the system are calculated. The optimization is re- 
garded as converged when the largest force and stress 
are smaller than 5 x 10~^ Ht/ae and 5 x 10"'' Ht/a^, re- 
spectively. We employ the maximally-localized Wannier 
function technique [ll|, which has been recognized to be 
a standard method for describing electronic orbitals in 
solid. There, Wannier functions are constructed from 
the Bloch wave functions belonging to the specified elec- 
tronic bands. In the present work we extract the transfer 
integrals between the localized Wannier functions con- 



structed for the soHd picene. 



III. RESULTS 

Pristine solid picene comprises molecular layers 
stacked in the c direction, as depicted in Fig. [TJ In 
each layer, picene molecules are in a herringbone arrange- 
ment, which means a unit cell contains two molecules. 
The crystal lattice is monoclinic with space group P2i. 
Its lattice parameters are reported to be a = 8.480, h = 
6.154, c = 13.515 in A, and the angle p = 90.46° be- 
tween a and c axes from a single-crystal X-ray diffraction 
experiment. [25I The samples used by Mitsuhashi et al. [l| 
had the same lattice parameters. 

Let us first classify the candidate positions into which 
dopants can be accommodated. A conventional organic 
insulator, solid pentacene, where each molecule has the 
same chemical formula C22H14 as picene but five ben- 
zene rings are connected on a straight line with zigzag 
edges, is known to accommodate dopants in its inter- 
layer region when doped with alkali elements, [26] with its 
c axis significantly expanded. By contrast, the c axis of 
picene is reported to shrink when potassium-doped with 
X = 2.9. [1] This suggests that the dopants have intruded 
into the intralayer region of solid picene. To determine 
whether intralayer insertion of potassium atoms is ener- 
getically more favorable than interlayer insertion, we per- 
formed in our previous work^ a structural optimization 
by locating a potassium atom on different positions on 
the ab plane of the pristine picene as initial geometries, 
and found that the intralayer insertion is in fact energet- 
ically favorable. In the present work, we hence explore 
optimized geometries of the dopants with intralayer in- 
sertion. 

The optimized lattice parameters and the internal co- 
ordinates of potassium atoms are summarized in Table HI 
The crystal symmetry of pristine picene was found to be 
preserved in all the optimized doped structures, despite 
significant deformations of each molecule from the planar 
shape along with significant rearrangements of molecules 
upon doping. The preserved symmetry contrasts with 
the situation in Kscoronene, for which we have shown 
that doping lowers the svmmetrv. p7| The geometries, 
the electronic band structures, the densities of states, 
and the Fermi surfaces for the optimized structures are 
shown for Kipicene in Fig. [21 for K2picene in Fig. [3l for 
Kspicene in Fig. HI and for K4picene in Fig. [5l In each 
case, the conduction band group comprises four bands, 
which have main characters of LUMO and LUMO+1 of 
an isolated picene molecule. Since the herringbone struc- 
ture has two molecules per unit cell, we have constructed 
four maximally-localized wave functions (WFs) from the 
four bands around the Fermi level of the individual struc- 
tures. For each of them, two WFs, wi^ and localized 
at each molecule with higher and lower orbital energies 
were obtained as shown in Fig. [6l 
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FIG. 1: (Color online) (a) Crystal structure of pristine solid 
picene (drawn with VESTA 24]), which comprises molecular 
layers stacked in the c direction. Solid lines delineate unit 
cells, (b) Typical first Brillouin zone of the doped picene op- 
timized in the present work. Due to the crystal symmetry, we 
have a* ± 6* and 6* ^ c*, while c* and a* are not orthogonal 
to each other. 



TABLE I: Optimized lattice constants a,b,c (in A), angle /3 
between a and c axes, volume V of a unit cell, and internal 
coordinates of dopants in parentheses for Kipicene with x = 
1-4. 



V 



Ki(A) 7.267 7.387 12.698 95.79 688.876 
(0.1820,0.3202,0.3615) 

Ki(B) 8.234 6.554 12.807 94.65 678.194 
(0.0101,0.3733,0.2382) 

K2 7.237 7.506 12.624 94.20 683.868 
(0.3916,0.3155,0.7333) 
(0.1829,0.3196,0.3935) 

K2K1 8.776 6.394 13.346 94.03 747.069 
(0.3243, 0.2864, 0.6301) 
(0.1806,0.2848,0.3225) 
(0.1941,0.2548,0.0331) 

K3(A) 7.421 7.213 14.028 104.53 726.848 
(0.3357, 0.2790, 0.8241) 
(0.2314, 0.2978, 0.5063) 
(0.1167,0.2792,0.1818) 

K3(B) 7.408 7.223 14.116 105.93 726.328 
(0.4355,0.2812,0.8178) 
(0.2247, 0.3007, 0.4932) 
(0.0105,0.2810,0.1755) 

K4 7.511 7.058 14.230 102.96 735.230 
(0.3039, 0.2585, 0.8344) 
(0.2635, 0.3334, 0.5928) 
(0.2441,0.3202,0.3477) 
(0.1377,0.2586,0.1147) 



A. Kipicene 

Let us start with Kipicene. Two structures, Ki(A) and 
Ki(B), have been obtained. Ki(A) appears when we put 
one potassium atom per molecule in the intralayer region 
as the initial position of the optimization, while Ki(B) 
appears when we put the dopant in inter layer. In both 
structures, the dopants end up with intralayer positions 
after the optimization. The insertion of the dopants into 
the intralayer region supports the result of our previous 
work. P!q| In Ki(B), with a total energy higher than in 
Ki(A) by 0.316 eV per molecule, the picene molecules 
are more strongly twisted than in Ki(A). Both systems 
have metallic electronic band structures. The width of 
the band group consisting of four around their Fermi level 
is 0.36 for Ki(A), and 0.52 eV for Ki(B). 



B. K2picene 

K2picene structure was obtained by putting one potas- 
sium atom in the intralayer region and another in the 



interlayer region as the initial positions of the optimiza- 
tion, where we confirmed that an additional atom doped 
into Ki picene again prefers the intralayer region by mov- 
ing to an intralayer position after the optimization. The 
obtained K2picene is an insulator, with a band gap of 
0.038 eV at X. The widths of the highest occupied and 
the lowest unoccupied band groups are 0.21 and 0.29 eV, 
respectively. By analyzing the projected density of states 
for the WFs, we find that each of the two band groups 
has comparable weights for the two types of the Wannier 
functions and wi. 



C. Kspicene 

For Kspicene we first obtaine Kspicene (A) structure by 
putting three potassium atoms in the intralayer region 
as the initial positions for the optimization. Andres et 
al [HI obtained the monoclinic Kspicene structure with 
a = 7.359,6 = 7.361, c = 14.018 in A and /3 = 105.71°, 
in which the potassium atoms are all in the intralayer re- 
gion. We have also tried a structural optimization start- 
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FIG. 2: (Color online) Geometries (large spheres: K atoms), electronic band structures, densities of states, and Fermi surfaces 
of the Ki structures. The origins of energy are set to the respective Fermi level. 
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FIG. 3: (Color online) Geometry (large spheres: K atoms), 
electronic band structure, and density of states of the 
K2picene. The origin of energy is set to the Fermi level. 



ing from their lattice parameters and atomic coordinates. 
Then the geometry resulted in another structure, which 
we call K3picene(B) here. Although the lattice param- 
eters of K3(A) and K3(B) are close to each other with 
the difference in the total energy being only 1 meV per 
molecule, they are distinct in the positions of dopants, as 
seen in Table |T] and in Fig. IH K3picene(B) has the lat- 
tice parameters and the dopant positions close to those 
of the K3picene structure obtained by Andres et al. [lH, 
and we therefore consider the two structures to be the 
same. K3(A) and K3(B) have metallic electronic band 
structures with the same width 0.6 eV of the band group 
consisting of four around their Fermi level. 

Since the superconductivity in potassium-doped picene 
was observed for x 3,[l|, [2| more elaborate explo- 
ration of geometries is required than for other concentra- 
tions. We have in fact found another structure, which 
we call K2Kipicene, by putting two potassium atoms 
in the intralayer region and one in the interlayer re- 
gion as the initial positions for the optimization. This 
structure is higher in the total energy than the K3(A) 
structure by 0.465 eV per molecule. However, the lat- 
tice parameters are closer to the experimental results [l| 
of a = 8.707, h = 5.912, c = 12.97 in A, and p = 92.77° 
for K2.9picene. The electronic band structure is again 
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metallic for the K2K1 structure. The four bands around 
the Fermi level, whose width is 0.51 eV, are close to the 
higher bands as compared with those in the K3(A) and 
K3(B) structures. 

Each of K3(A) and K3(B) has multiple Fermi sur- 
faces comprising a pocket and complicated, connected 
surfaces. The existence of a Fermi pocket implies three- 
dimensional electron transfer integrals despite the layered 
structure. On the other hand, K2K1 has one-dimensional 
(a pair of planar) Fermi surfaces along with connected 
surfaces. We can also note that, while the Fermi sur- 
face of K3Coronene[23 also consists of multiple sur- 
faces, the topology of the Fermi surfaces for K3Coronene 
(consisting of a pair of quasi-one-dimensional, planar 
surfaces and another pocket) differs from those found 
here for K3picene(A) and K3picene(B) (consisting of 
multiple surfaces of two- or higher dimensional charac- 
ters), and for K2Kipicene (consisting of a pair of quasi- 
one-dimensional, planar surfaces and another, multiply- 
connected one). 

The largest interlayer transfer integral of the WFs in 
K3(A) (K2K1) is 13.6 (19.6) meV, while that in undoped 
picene is 14.6 meV. The larger interlayer transfer inte- 
grals in K2K1 than in K3(A) originate in the presence 
of dopants in the interlayer region. The widths of the 
bands around the Fermi level in these doped structures 
are larger than in undoped picene (0.39 eV[Io[), though 
K3(A) has the smaller interlayer transfer integrals. It 
implies that their larger band widths come from the en- 
hanced intralayer transfers of electrons due to the hy- 
bridization of the wave functions of potassium atoms and 
the molecular orbit als. In fact, the largest intralayer 
transfer integral in K3(A) (K2K1) is 82.7 (70.6) meV, 
which is much larger than that in undoped picene, 57.9 
meV. The intermolecular transfer integrals between the 
WFs of K3(A)picene and K2Kipicene are shown in Fig. 
[71 It is seen that K3(A) allows for the transfers between 
more distant WFs in the a direction, which leads to its 
more two-dimensional and wider band dispersion than 

K2K1. 



D. K4picene 

We have obtained the K4picene structure by putting 
three potassium atoms in the intralayer region and one 
in the interlayer region as the initial positions of the 
optimization, because the three dopants sitting on ev- 
ery other benzene rings out of five rings in each picene 
molecule may tend to push the fourth dopant to the in- 
terlayer region. In the optimized structure, all the four 
dopants end up with sitting in the intralayer region, but 
two of them are rather close to the interlayer region due 
to the crowded space. The system is an insulator with a 
band gap of 0.36 eV. The highest occupied band group 
consists of two bands, separated from the second highest 
group of two bands by 0.01 eV. The widths of the higher 
and lower band groups are 0.17 and 0.20 eV, respectively. 



By analyzing the projected density of states for the WFs, 
we found that the higher and lower occupied band groups 
consist mainly of and respectively. This contrasts 
with the situation in another insulator, K2picene. 



E. Discussions 

For each concentration x, the lowest-energy structure 
shown above contains the dopants in the intralayer re- 
gion. This behavior is consistent with the observed 
tendency [2] of the c axis of K^^picene, in which the c 
axis shrinks monotonically with increasing x, implying 
that the dopants move into the intralayer region regard- 
less of the concentration. This is to be contrasted with 
Rba^picene, whose c axis was observed to expand up to 
X = 2 and shrink afterwards. 

In Fig. [6] we have displayed the maximally-localized 
Wannier functions (w\^ and w\ with higher and lower 
orbital energies, respectively) in the optimized struc- 
tures for all the cases of Kipicene to K4picene. We 
can see that the WFs have basically the same features 
for all the structures in that all the WFs have LUMO 
and LUMO+1 orbital characters, although small but sig- 
nificant amplitudes found around the potassium atoms 
whose magnitudes depend on the doping level x. More 
precisely, we can show that and w\ have charac- 
ters of the linear combinations of LUMO and LUMO+1: 
^h,i = (/>LUMO ± 0LUMO+1, as in the undoped picene p^. 

However, a key message in the present band struc- 
ture results is that, although the orbital characters of 
the WFs are common to all the optimized geometries, 
their band structures very much differ. This is not sur- 
prising, since the molecular shape as well as the herring- 
bone arrangement of molecules distort significantly with 
doping. More importantly, the lobes of molecular WFs 
extend outward the molecule via the potassium atoms, 
with significant hybridizations of WFs around the potas- 
sium atoms. Specifically, we find that the small contribu- 
tions from the potassium wave functions have a p orbital 
character, as seen in Fig. [6l 

In other words, the rigid-band picture is invalidated 
for multifold reasons, i.e., through the rearrangement 
and distortion of molecules, but also via hybridization 
of the potassium p orbitals with the molecular orbitals 
of picene. Indeed, for all the optimized geometries we 
have also calculated their band structures (not shown) 
with potassium atoms removed but the doping level is 
fixed by shifting the chemical potential to find that the 
band structures are drastically changed from those of the 
doped systems. As a consequence the resultant Fermi 
surface exhibits a variety of multiband structures which 
take diverse topology for Kipicene and K3picene. This 
obviously provides an intriguing starting point for explor- 
ing mechanisms of superconductivity as a future work. 

Final point of interest, about the molecular species de- 
pendence, is the following: we can note that the cal- 
culated band structures for all x (including x = Op^of) 
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FIG. 4: (Color online) Geometries (large spheres: K atoms), electronic band structures, densities of states, and Fermi surfaces 
of the Kspicene. The origins of energy are set to the respective Fermi level. 



have significant dispersion along the c direction as well 
as along the a and b directions, which indicates isotropic 
components in the electron transfers. If we compare this 
with the situation in K^coronene, we find that the band 



dispersion is more anisotropic in both undoped coronene 
and K 3 coronene [13, which reflects the differences be- 
tween the intermolecular distances in each direction. 
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FIG. 5: (Color online) Geometry (large spheres: K atoms), 
electronic band structure, density of states, and Fermi sur- 
faces of the K4picene. The origin of energy is set to the Fermi 
level. 



IV. SUMMARY 

We have systematically performed a full structural 
optimization for K^^picene, covering the doping level 
X = 1-4. The crystal symmetry of pristine picene of 
P2i was found to be preserved in all the optimized 
structures despite the deformation and rearrangements 
of the molecules upon doping. The optimized systems 
for Kipicene and Kspicene have metallic band struc- 
tures, while band gaps open around their Fermi levels 
for K2picene and K4picene. In each of the optimized sys- 
tems, the electronic bands in the vicinity of the Fermi 
level are derived from LUMO and LUMO+1 of the aro- 
matic molecule as in the undoped picene. We found that 
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the presence of dopants affects the electronic properties 
of all the doped structures. The electronic structures 
obtained with the local density approximation for the 
optimized structures reveal that the rigid-band approxi- 
mation is invalidated due to the rearrangement and de- 
formation of molecules and the accompanied charge re- 
distribution between the molecules and dopants. The 
resultant Fermi surface exhibits a variety of multiband 
structures which take diverse topology for Ki picene and 
Kspicene. 

In all the optimized systems, the widths of the LUMO- 
and LUMO+l-derived bands do not exceed 0.6 eV, which 
are much smaller than the experimental electronic corre- 
lation energy U ~ 0.85 eV,[i5[ indicating that the mate- 
rial may be a strongly correlated electronic system. Fur- 
ther investigations for clarifying the material properties 
and the mechanism of superconductivity in doped picene 
should thus take into account the magnetic instability 
and the strong Coulomb interaction, as pointed out by 
Giovanetti and Capone[13] and Kim et al. Q 

Rbs.i picene and Cai.spicene have also been reported 
to become superconducting with ^ 7 K,[2| which indi- 
cates that the transfer of three electrons from the alkali 
atoms to a picene molecule can lead to superconductivity 
in alkali-doped picene. Theoretical exploration of their 
possible geometries is needed, where Rba^picene is ex- 
pected to accommodate dopants in the interlayer region 
due to a larger atomic radius of Rb than in K^^ picene. 
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